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We consider the analogy between the topological phase transition which occurs as a function of spatial coordinate 
on a surface of a non-trivial insulator, and the one which occurs in the bulk due to the change of internal parameters 
(such as crystal field and spin-orbit coupling). In both cases the system exhibits a Dirac cone, which is the universal 
manifestation of topological phase transition, independently on the type of driving parameters. In particular, this leads 
to a simple way of determining the topological class based solely on the bulk information even for the systems with 
translational symmetry broken by atomic disorder or by strong electron correlations. Here we demonstrate this on ex- 
ample of the zinc-blende related semiconductors by means of the ab-initio fully-relativistic band structure calculations 
involving the coherent potential approximation (CPA) technique. 
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1 Introduction. The remarkable ability of the semi- 
conductors with topologically non-trivial band structure 
(for the overveiw see e. g. Ref. (H) to exhibit the non- 
dissipative quantized surface spin-current attracts great at- 
tention both in theoretical physics, material science and 
spintronics engineering. This spin-current is represented 
by a pair of opposite charge currents carrying opposite 
spin states. Each charge current is induced by the intrin- 
sic spin-orbit interaction, which can be viewed as spin- 
dependent effective magnetic field. These effective mag- 
netic fields quantize the electron orbits into the Landau 
levels, which leads to a quantization of the resulting spin- 
current in units of (1/2 + n) (e 2 /h) . The absence of dissi- 
pation holds as long as the structure of the Landau levels 
is preserved which is guaranteed by the time-reversal sym- 
metry. In contrast to the bulk, the non-compensated spin- 
current can be observed and detected on the surface, due to 
the space-reversal symmetry break while keeping the time- 
reversal symmetry. On the surface the system exhibits the 
Dirac cone, which can be observed directly by the angular- 
resolved photoemission technique. This property of topo- 
logical insulators is typically exploited in the band struc- 
ture calculations in order to probe the non-trivial nature of 
the material. On the other hand, the reliable surface cal- 
culations most often appear to be extremely resource de- 
manding. Alternatively, since the topological edge states 



emerge from the non-trivial bulk band structure, it must 
be a way for determining the topological class based on 
pure bulk information. Indeed, it is possible to calculate 
the topological invariant Z2 based on the parity analysis 
of the occupied eigenstates |2,3]. Obviously, this proce- 
dure cannot be applied in situations when the eigenstates 
do not exist due to the absence of the translational sym- 
metry. The reasons for such symmetry break are typically 
the non-periodic boundary conditions, atomic (chemical 
or structural) disorder, or the intrinsic electronic disorder 
caused by strong dynamical correlations. The hint is that 
in the first case, i. e. on the surface of topological insula- 
tor one always observes the Dirac cone. At this moment 
it is important to understand that the surface is not just a 
complicated geometrical object with reduced dimension- 
ality, potential step etc., but it is a transition between dif- 
ferent bulk phases. The parameter of this transition is sim- 
ply the spatial coordinate. In this context the Dirac cone 
is not an ingenuous surface phenomenon, but is rather a 
fundamental manifestation of the transition between two 
distinct topological phases in certain parameters space. In- 
deed, the Dirac cone can occur in the bulk band structure as 
well: all materials which appear at the borderline between 
topologically non-trivial and trivial phases exhibit the bulk 
Dirac cone (H . This fact indicates that for materials with 
the translational symmetry broken by disorder the topolog- 
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Figure 1 (color online) Left: crystal structure of the zinc- 
blende binary YZ compound (space group No. 216). Yel- 
low an red spheres represent the transition-metal Y and 
main-group elements Z. Right: by adding one more tran- 
sition element X (blue sphere) one obtains ternary XYZ 
Heusler compound. 



ically non-trivial state still can be determined based on the 
bulk information. 

In the present work we would like to illustrate this fact 
by first-principles calculations taking into account the ran- 
dom disorder by means of the Coherent Potential Approxi- 
mation (CPA) alloy theory . Simultaneously, the CPA 
serves as adiabatic transition technique, which allows us to 
model the smooth transformation of one material into an 
other and observe the details of the topological phase tran- 
sition in the bulk. This transformation has also a physical 
relevance, i. e. it directly describes the effect of the ran- 
dom chemical substitution. As a real-case example we will 
consider the class of zinc-blende related structures, which 
contain a great manifold of non-trivial and trivial semicon- 
ductors in which the details of their band structures are 
well-known. 

2 Topologically non-trivial systems with zinc- 
blende structure. Topological semimetals (or gapless 
semiconductors with topologically non-trivial band struc- 
ture) represent a rich and at the same time rather simple 
class of materials, which allows to demonstrate the general 
trends specific for the class of the 2D topological insu- 
lators as a whole. The simplest systems in this class are 
the binaries with zinc-blende structure (see Fig. [U as e. g. 
HgTe [7,8]. It is relatively easy to search for the nonmag- 
netic semiconductors in this class by picking the pairs of 
elements with average number of valence electrons per 
atom n va i = 6 |4, 9 ]. Without restriction one can add here 
the whole class of ternary Heusler semiconducting materi- 
als (so-called half-Heusler semiconductors) which possess 
the same space group No. 216 and are different from bina- 
ries by having one more transition metal element l4ll9l[T0l . 
As it follows from these studies, all topologically-related 
physics of zinc-blende like compounds is contained in the 
vicinity of the T-point: all trivial semiconductors of this 
type exhibit a direct non-zero band gap, whereas all non- 
trivial materials are gapless. It is surprising, but there was 
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Figure 2 (color online) Relativistic band structures calcu- 
lated for the ordered YZ binary zinc-blende (ZnS, CdTe, 
TIP, HgTe and CdPo) and related ternary XYZ Heusler 
(LuPtSb) materials. Red, green and blue colors refer to 
the orbital symmetries sy 2 (carried by the transition ele- 
ment Y),py 2 and p 3 / 2 (carried by the main group element 
Z), respectively. The gray arrows indicate the transition 
from trivial semiconductors ZnS, CdTe towards topolog- 
ical phase transition marked by the Dirac cone (LuPtSb), 
and further into the class of the non-trivial gapless semi- 
conductors TIP, HgTe and CdPo. Calculations are per- 
formed by the PY-LMTO method El. As the LDA El 
form of the exchange-correlation potential was used, the 
band gaps of trivial semiconductors (ZnS, CdTe) are un- 
derestimated. 



no exception found yet. Here we also illustrate the differ- 
ence between topological and trivial materials within this 
class in terms of the band structure, by giving few examples 
shown in Fig. [2] The trivial band structure is characterized 
by the empty sy 2 orbital of the transition metal atom Y and 
the fully filled p-shell of the main-group element Z, split 
into p 3 / 2 and py 2 by spin-orbit coupling. These compounds 
exhibit a non-zero direct band-gap at the Fermi energy and 
are represented by the light-element semiconductors with 
average nuclear number (Z) roughly less than 60: e. g. ZnS 
((Z) = 23), ZnSe ((Z) = 32), CdTe, ScPtSb ((Z) = 50); 
with certain exceptions, as e.g. LuNiBi ((Z) w 60.7) E). 
In the non-trivial systems the p 3 / 2 orbital is half-filled, for 
this reason it sticks to the Fermi energy and within cubic 
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symmetry the system becomes zero-gap semiconductor. 
At the same time the py 2 and sy 2 orbitals are fully occu- 
pied. This so-called band inversion changes the parity of 
the wave function which is directly related to the topolog- 
ical Z2 invariant. Such materials are typically heavy with 
(Z) > 60: ZnPo ((Z) = 57), HgTe and CdPo ((Z) = 66), 
YPtBi ((Z) w 66.7), LuPtBi ((Z) « 77.3) Hill; with cer- 
tain exceptions, as e. g. was found in present calculations 
for TIP ((Z) = 48). Interesting case is represented by the 
borderline compounds, such as YPdBi, YPtSb ((Z) = 56), 
or LuPtSb ((Z) w 66.7) 1 10], which exhibit the Dirac cone 
in the bulk, i.e. can be easily switched between topological 
or trivial states by small variations of the stoichiometry or 
the lattice constant (H. 

Obviously, there is also a continuum of alloys, which 
can be obtained by combining various ordered zinc- 
blende (Yi- y Yy)(Zi- z Z' z ) or half-Heusler materials 
(X 1 - x X , x )(Y 1 - y Y , y )(Z 1 - z Z , z ). By varying the ratios of 
elements allows to tune the properties of the system con- 
tinuously in a wide range. For example, by substituting 
Hg with Cd in (Hgi_ cc Cd ;r )Te series one obtains the zero- 
gap materials for x < 0.17 and non-zero band gap for 
x > 0.83 (HI, which correspond to the non-trivial and 
trivial regimes, respectively. 

3 Chemical disorder in terms of the CPA alloy 
theory. In the context of the first-principles description 
of the electronic structure, the most practical way to treat 
the random disorder is by using the so-called Coherent Po- 
tential Approximation (CPA) 0151. In fact, until now the 
CPA (and its extensions) remains the only technique which 
incorporates the effects of the energy-dependent shift and 
lifetime broadening - the essential features of the electron 
localization caused by chemical disorder, which are not ac- 
cessible within other theories, as e. g. the VCA [ITU or su- 
percell calculations. 

The idea of CPA can be sketched by considering the 
example of binary alloy A^Bi-^ (0 < x < 1), with each 
lattice site occupied either by atomic sorts A or B with 
probabilities x or 1 — x, respectively. For simplicity we let 
the atomic potential to be characterized by a single energy 
level, 6a and ee, for each sort respectively. In the second 
quantization form the Hamiltonian of A^Bi-^ alloy Hab 
can be written as 

Hab = ^2 eka k a k + ^2 e i a l a i • (!) 

k i 

The first term represents the sum over the valence states 
indexed by their wave vectors k with kinetic energies e&. 
The second term is the sum over the lattice sites i; ei is 
the randomly fluctuating variable, equal to ca or cb with 
probabilities x or 1 — x, respectively. 

The straightforward solution for such Hamiltonian rep- 
resents a complicated many-body problem. To make it 
practically accessible, one needs to find an adequate mean- 
field approximation, i. e. to substitute the real disordered 
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Figure 3 (color online) Mean-field treatment of random 
disorder. Random binary A^Bi-^ alloy (A and B atoms are 
shown as red and blue spheres) is substituted by the regular 
array of the single effective atomic sort C (purple spheres). 
The effective potential ec is adjusted in order to exhibit 
the same properties as the original disordered binary: for 
example, the same density of states. 

system A^Bi-^ by the effective regular system C 

H A b ~ H c = ^ e k a\a k + ^ e c a - a* 

k i 
k i 

which emulates the same properties. This idea is sketched 
on Fig. [3] The quality of such emulation is a matter of how 
the effective potential ec is adjusted. It is easy to show, 
for example, that the straightforward averaging over the 
single-site potentials misses much of important physics. In- 
deed, since both ca and cb are the real-valued quantities, 
the effective regular potential ec = x • ca + (1 — x) • cb 
will be real- valued as well. Thus, the solutions of the cor- 
responding Schrodinger equation will be the delocalized 
Bloch waves which propagate without scattering and al- 
ways lead to an infinite residual conductivity in case of a 
metal. 

The CPA is based on a more rational choice of the 
mean-field variable by using the one-particle Green's func- 
tion Gab • Since it represents a response of the whole sys- 
tem to the (5-like perturbation, i. e. 

[e - Hab(t)} G AB (r, r', e) = 5(r - r') , (3) 

it can be seen as a weighted superposition of responses 
from all configurations Q present in the system: 

Gab = ^2x q Gq , y^x Q = l. (4) 

Each configuration Q can be presented as a perturbation of 
the effective regular system C, i. e. it can be defined using 
the Dyson equation: 

Gq = Gq + Gq ^2 ( €i ' n ~ ec ) 

i 

= Gq + GqTqGq , (5) 
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Figure 4 (color online) Single-site CPA. Random binary 
Ai-^B^ alloy (A and B atoms are shown as red and blue 
spheres) is approximated as a superposition of two random 
impurity Green's functions (the impurities of types A and 
B are embedded in the regular host of the effective atoms 
C shown in purple) weighted with corresponding concen- 
trations x and 1 — x. 



where is equal ca or cb, depending on the configura- 
tion Q\ the perturbation with respect to the regular system 
C is fully described by the term 

= Ej - eg) 

Thus, the requirement 

Gab = Gc 



(6) 



(7) 

is equivalent to the cancellation of all partial perturbations 
in a total sum 



(8) 



Q 



which completes the self-consistent set of equations, 
needed to determine eo Since the definition © involves 
the Green's function Gc, which is complex- valued and 
energy-dependent, the effective potential ec appears to be 
complex- valued and energy-dependent as well. Its real part 
Re ec defines the energy position of the electronic state, 
whereas the imaginary part Imec is responsible for its 
finite lifetime broadening caused by disorder. The latter 
is the necessary ingredient for the calculation of residual 
conductivity (except of the Berry phase contribution fT5l 
which is a pure band structure effect). At the same time, 
the effective system formally remains to be translational 
invariant, thus once ec is determined, we can restrict all 
further calculations to the primitive unit cell. 

In order to simplify the set of the mean-field equations 
(E [6l [8]) we assume further that the only important con- 
figurations are those where the atoms are placed fully ran- 
domly, i. e. without certain nearest environment preference. 
In this case, only two configurations are left: random sin- 
gle impurities of type A and of type B, embedded in the 
C host (see Figure©. Thus, the shortcoming of the single- 
site CPA is the absence of local environment effects. The 
latter can be, however, taken into account by its non-local 
extensions fT6l[T7,181. On the other hand, for the case of 



random isovalent substitution, the single-site CPA remains 
a quite reasonable approximation even in the diluted limit, 
since the isovalent atoms intermix without extra environ- 
mental preference. The only condition is to preserve the 
zinc-blende structure which is again guaranteed by isova- 
lency. 

In the following all electronic structure calculations 
were carried out using the fully-relativistic Korringa- 
Kohn-Rostoker (KKR) Green's function method imple- 
mented within the SPR-KKR package 1 19 ] employing the 
density functional theory framework. Due to the multiple- 
scattering construction of the Green's function l20l , the 
method provides a suitable base for mean-field approaches 
like CPA. Since the electronic structure of disordered 
materials in general does not show well defined Bloch 
eigenstates, their electronic structure is described by the 
Bloch- spectral function (BSF) defined by the imaginary 
part of the alloy Green's function that is diagonal in mo- 
mentum space (20]. The exchange and correlation was 
treated using the Vosko-Wilk-Nusair form of the local 
density approximation (LDA) fT2l . LDA is known for its 
typical underestimation of the band gaps of semiconduc- 
tors. However, in case of our study this issue is not critical, 
since the observed trends are qualitatively well described 
even with LDA. 

4 Topological bulk phase transition induced by 
chemical substitution. In the following we will track 
the evolution of the BSF driven by random chemical sub- 
stitution along three distinct adiabatic paths, as depicted in 
Fig. first by going between the pair of trivial, second - 
between the pair of non-trivial materials. In the third case 
we span the path from the non-trivial to the trivial one in or- 
der to study the details of the topological phase transition. 
In order to keep the conditions of adiabatic transition we 
perform the random substitution in isovalent manner, by 
keeping the chemical bonding rules and the crystal struc- 
ture: i. e. YZ (Yi-sYj.XZi-sZj.) Y'Z'. For the alloy 
series the lattice constant is approximated linearly between 
the experimental constants for the pure materials. 

As an example for the first case, we select two trivial 
materials, namely CdTe and ZnSe which are well-known 
semiconductors with non-zero band gaps. Both systems ex- 
hibit the same band structure features (Fig. [6]): the conduc- 
tion band is formed by the sy 2 - symmetry (5s for Cd and 4s 
for Zn), whereas the top of the valence band carries the p 3 / 2 
symmetry (bp and 4p orbitals for Te and Se, respectively). 
As it follows from the BSFs for the intermediate composi- 
tions, the random substitution of Cd by isovalent Zn and Te 
by isovalent Se does not change the energy distribution of 
the orbital characters, thus all the alloys in between exhibit 
similar band structures. On the other hand, one can easily 
notice their principle difference from the ordered materi- 
als: for < x < 1 the symmetry analysis of eigenvalues 
cannot be applied directly, since there are no pure Bloch 
eigenstates, i. e. in the disordered regime the (5-like poles 
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Figure 5 (color online) Schematic illustration of the topo- 
logical phase transition via substitution. Light red and blue 
areas distinguish between manifolds of topologically non- 
trivial and trivial materials, respectively. The correspond- 
ing band structures are schematically given in the vicin- 
ity of the red and blue spheres which indicate certain 
compounds from topological and trivial manifolds, respec- 
tively. The arrows indicate the transformation of one mate- 
rial into the other, e. g. by chemical substitution or by over- 
going through the interface. The Dirac cone indicated by 
light purple lines always appears for the borderline case. 



of the BSF turn into overlapping Lorentzians, spread over 
the whole /c-space. Despite that we still can trace the evo- 
lution of their maxima in analogy to dispersion relation for 
the pure material, the only strict fact which tells us that all 
these materials belong to the same topological class is the 
absence of the Dirac-like dispersion within these series. 

In the second case we transform the topologically non- 
trivial HgTe into another non-trivial zinc-blende system 
CdPo. The Po zinc-blende compounds (CdPo, ZnPo) were 
not considered earlier in the context of topology, how- 
ever their crystal and band structures are well-known long 
ago H2 1112211 . As the presently calculated band structures 
show, CdPo is a gapless semiconductor with the bands 
alignment equivalent to HgTe (see Fig. 0. As it follows 
from the calculated series shown on Fig. [71 the situation is 
in principle very similar to the first case, with only dif- 
ference that all materials belong to the non-trivial topo- 
logical class, i. e. (Hg^Cdi-^XTe^Poi-^), (0 < x < 1) re- 
main gapless semiconductors and none of them exhibits 
the Dirac cone. On the other hand, it is interesting to ad- 
mit, that due to different energy placements of the s-state 
in pure HgTe and CdPo compounds one observes the cone- 
like feature for the composition x « 0.3 at about -1.5 eV. 
Indeed, in case of HgTe the s-state is situated below both 
pi/ 2 and p 3 / 2 and in case of CdPo - right in between (see 
Fig. 0. By substituting Te to Po, the effective spin-orbit 
split increases by pushing the py 2 -character down in en- 
ergy until it interchanges its position with s-character state 
at x « 0.3. Since this re-ordering of the states is not con- 
nected with their re-occupation, the parity of the total wave 
function does not change and the rest of the alloy series re- 
mains in the same topological class. 



In the last case the adiabatic path is going from 
the non-trivial CdPo to the trivial ZnSe semiconduc- 
tor. As it follows from BSFs calculated for the series 
(Cdi-^Zn^XPoi-^Se^) (Fig. [8]), we indeed encounter cer- 
tain composition (x c ~ 0.72) characterized by the bulk 
Dirac cone. All other alloys with x < 0.72 are gapless 
semiconductors and topologically non-trivial, whereas 
those for x > 0.72 are in the trivial regime, with the non- 
zero band gap. We briefly describe the most important 
steps of the BSF evolution in the vicinity of the r point: 
the Po 6p-band is split due to spin-orbit coupling into a 
p 3 / 2 -state (of i~8 symmetry in the pure case) sitting right at 
the Fermi energy and a py 2 -state (of I7 symmetry) which 
is about 3 eV lower. The Cd 5si/ 2 band (of i~6 symmetry) 
is placed by the crystal field in between, at about -2 eV. 
Since by substitution we effectively decrease the spin-orbit 
coupling and increase the crystal field, both pu 2 and su 2 - 
like states rise in energy by shifting towards the Fermi 
level. At the "critical" composition x c = 0.72 the system 
undergoes a topological phase transition manifested by the 
Dirac cone formed by a mixture of p- and s-symmetric 
states. 

Due to the extremely deep energy position of the Cd 
5-shell in CdPo compound (approximately at -2 eV), the 
topological phase transition occurs within rather ZnSe-rich 
regime x c = 0.72 > 0.5, which is not so severely disor- 
dered as in case of x = 0.5. Indeed, as it follows from 
the similar adiabatic path connecting non-trivial HgTe with 
trivial ZnSe semiconductor (see Fig. [9]), the borderline con- 
centration is shifted towards non-trivial regime since the 
s-symmetric state of the pure non-trivial alloy component 
is situated closer to the Fermi energy. Despite that the 
borderline composition is then shifted closer to the max- 
imally disordered system x c = 0.63 « 0.5, the Dirac cone 
exhibits rather well-defined B loch-like states compared to 
the rest of the band structure. 

Obviously, that the edge states with linear dispersion at 
the surface of a topological insulator and the states form- 
ing the Dirac cone in the bulk manifest the same transition 
phenomenon. The qualitative difference is only that due 
to the remaining effective translational symmetry (the dif- 
ference is that the electrons are localized by disorder and 
are not Bloch-like), we observe two replica of Dirac cones 
with opposite spins on top of one another, whereas due to 
the break of space-reversal symmetry at the surface only 
a single cone remains, which exposes the adiabatic spin- 
current. 
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topologically non-trivial and characterized by zero-gap, whereas x > 0.63 are trivial and exhibit a real band gap, which 
increases towards x = 1. The composition at x = 0.63 is characterized by the clear Dirac cone centered at the Fermi 
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